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METHOD FOR CHARACTERISING AN OPTICAL FIBRE LINK 

10 Field of the invention 

[0001] The present invention is related to a method 

for characterising an optical fibre link by analysis of a 
Rayleigh backscatter signal, from which the spatial 
distribution of the beat length, the coupling length and 

15 the PMD (Polarization Mode Dispersion) along the link can 
be determined. 

State of the art 

[0002] The use of dispersion-shifted and dispersion 

20 compensating fibres has minimised the effect of chromatic 
dispersion on the bandwidth of an optical link. 
Polarisation mode dispersion (PMD) has therefore become the 
most serious limiting factor in high-speed optical 

♦ 

communication systems. There exist several measurement 
25 techniques of PMD, but they only allow for measuring the 
global value and do not give information about the 
distribution of PMD along the fibre length. Measuring the 
spatial distribution of PMD is, however, important. Such a 
measurement would indeed allow to locate the bad trunks in 
30 an optical link, causing high PMD value. This can be 
essential in the frame of a network maintenance and 
upgrade . 

[0003] The PMD of a fibre depends on two parameters: 

the beat length (Lb) and the coupling length (L c ) . The beat 
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length depends on the birefringence of the fibre and the 
coupl length is related to the phenomenon of 

polarisation mode coupling. Measuring the spatial 
distribution of PMD therefore consists of measuring the 
5 d±stribution of both the beat and coupling lengths. 

[O0 04] Most of the techniques proposed for 

measurement of the distributed PMD are based on 
polarisation-optical time domain ref lectometry (POTDR) . The 
concept of POTDR was introduced by A. Rogers in the 80 's 

10 (see A . Rogers, s Polarization-optical time domain 
r&f lectometry: a technique for the measurement of field 
distributions' , Applied Optics, vol. 20, pp. 1060-1074, 
19 81) . It basically consists of measuring the polarisation 
properties of the Rayleigh backscatter signal when an 

15 optical pulse propagates down the fibre. Some POTDR-based 
techniques have since been developed for PMD 
characterisation. Most of them allow a measurement of the 
beat length distribution along an optical link and, 
therefore, do not allow the complete determination of the 

20 PMD . However, Sunnerud presented a POTDR-based set-up 
enabling the measurement of the accumulation of PMD along 
the fibre (see H. Sunnerud et al., x Polarization -Mode 
Dispersion measurements along installed optical fibers 
us±ng gated backscattered light and a polarimeter' , J*, of 

25 Lightwave Technol . , vol. 18, pp. 897-904, 2000). This 
method requires POTDR measurements at several wavelengths. 
A technique based on the measurement of the degree of 
polarisation by using a POTDR set-up has also been 
presented in ^Distributed detection of high-PMD sections on 

30 installed fibers using a polarization-OTDR' , M. Leblanc, 
Proc. OFMC'01, pp. 155-162, 2001, and is related to the 
patent application US 2003/0174312. This technique does not 
allow quantification of the PMD, but enables to obtain 
distributed information about the level of PMD. Moreover 
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this technique requires the complete measurement of the 
state of polarisation of the backscattered signal . More 
recently, Galtarossa described another technique for the 
PMD measurement for which the coupling length is determined 
5 from the correlation of the birefringence vector (see 
'Measurement of birefringence correlation length in long, 
single-mode fibers', Galtarossa. et al. , Opt. Lett., vol. 
26, pp. 962-964, 2001) , which also requires the measurement 
of the polarisation states of the backscattered signal. 

10 Moreover, this technique requires a quite small duration of 
the pulses launched in the optical fibre: 5-10 ns. 
[0005] In x Polarization mode dispersion mapping in 

optical fibers with a Polar ization-OTDR' , M. Wuilpart, G. 
Ravet, P. Megret and M. Blondel, Photon. Technol . Lett., 

15 vol. 14, pp. 1716-1718 , 2002, the authors describe an 
analysis of the POTDR signal which allows a mapping of the 
PMD along an optical fibre link by quantifying the PMD on 
each fibre of the link. The main advantage of that 
technique is that it does not require the complete 

20 measurement of the backscattered polarisation state 
evolution and only uses a linear polariser at the fibre 
input, which is quite simple to implement. The 
determination of the beat length, the coupling length and 
finally the PMD is based on the analysis of the statistical 

25 properties of the extrema present in a POTDR trace. 

[0006] In said paper by Wuilpart et al . the spatial 

distribution of the linear birefringence 8 is characterised 
by a Rayleigh statistical distribution along the fibre 
length and it is assumed that the rate of change of the 
30 birefringence angle q is described by a Gaussian 

distribution with zero mean and standard deviation a 
depending on the polarisation mode coupling strength. The 
circular birefringence is neglected. The corresponding 
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measurement set-up is shown in Fig.l. The OTDR pulses (10 
ns) externally modulate a 1553 nm DFB laser via a pulse 
generator and an electro-optic modulator (EOM) . After 
amplification, the pulses are launched into the fibre 
5 throixgh an acousto-optic modulator (AOM) , which suppresses 
the amplified spontaneous emission noise of the erbium 
doped, fibre amplifier between two successive pulses. A 
linear polariser is placed at the fibre input and a 
polarisation controller (CP) is used to obtain the maximum 
10 power after the polariser. 

[0007 3 Fig. 2 shows three simulated POTDR traces 

corresponding to the POTDR arrangement and obtained for 
different values of the beat length and an identical a 
value . One can easily observe that the number of maxima 

15 strongly depends on L B - 

[0008] The calculation of the beat length is based 

on the measurement of the number of maxima of the POTDR 
trace. Fig. 3, which has been obtained by simulation, 
presents the number n of maxima per unit of length in 

20 function of the beat length for different values of a. This 
figure shows that n varies not only with L B but also with 
the strength of the polarisation mode coupling (a) . In 
order to estimate L B from n, the curve corresponding to a = 
10 (thie dashed curve in Fig. 3) is used. This choice indeed 

25 minimi. ses the maximum possible error if one assumes that a 
is smaller than 15 degrees/m. Consequently, this method 
allows only to give an approximate value of the beat length 
Lb. A relationship between the coupling length L c and a can 
be derived and one can obtain: 

T - 0 879 

30 L C ~ — 

[0009] A statistical approach is then used for 

determining the coupling length from the POTDR trace. Fig. 
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4 shows different POTDR traces that were obtained by 
simulation for a beat length of 25 m and for different 

values of a: 2, 4 and 8 degrees/m, respectively. One can 
observe that the lower envelope of the curves varies in a 
5 smoother way when a decreases and therefore when L c 

increases. Defining APi as the absolute value of the 
difference in ordinate between two successive minima of the 
POTDR trace P one can write : 



10 where z± and z± + i correspond to the distances of the i and 
the i + l th minimum, respectively. 



where < . >i denotes the mean value for i varying from 1 to 
15 the number of minima of the POTDR trace. £ is a measure of 
the speed variation of the lower envelope of the 
backseat tered signal along the fibre length. Fig. 5 
(obtained by simulation) shows the evolution of £ according 
to a for different values of the beat length L B . The values 
20 are here expressed in dB because £ has been calculated 
directly from the POTDR trace. 

[0011] The procedure for the determination of L c is 

the following. From the POTDR trace, the number of maxima 
and each APi are measured. An approximate value of the beat 

25 length is then deduced by the method previously described 
and £ is calculated by averaging APi. Finally a can be 
deduced from £ by using the curve of Fig. 5 corresponding 
to the correct beat length and then L c is calculated by 
means of the relation L c =0 . 879/cr 2 stated above. Because Lb 

30 is an approximate value, the calculated coupling length L c 
will also be an approximation. 



A^=|P(^)-P(z m )| 



[0010] 
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[0012] From the measurement of the beat length and 

the coupling length, it is easy to determine an approximate 
value of the PMD value by 



PMD 1 =- 
2 



f X ^ 



2 



v L c 

K c J 



5 where c is the light velocity in a vacuum, X the 
measurement wavelength and L the fibre length. It should be 
emphasised that the final PMD value is approximate: this is 
due to the approximation done on the beat length 
measurement . This process can be repeated for each 
10 contribution in the POTDR signal corresponding to the 
different fibres of the optical link: it is therefore 
possible to map the PMD along this link. 

[0013] The method as described above has some 

important drawbacks and limitations. It only yields 

15 approximate values of L B , L c and the PMD. Further, the 
presented POTDR traces were assumed ideal. It supposes 
ideal characteristics of the various components of the 
experimental set-up. In practice, the set-up is 
characterised by some imperfections, which cause a 

20 distortion of the ideal POTDR signal. Therefore the method 
as such cannot be directly applied. The POTDR signal is 
also affected by noise: detector noise and residual 
coherence noise, which adds a series of minima and maxima 
to the POTDR trace. The final value of the measured PMD can 

25 therefore be erroneous. 



Aims of the invention 

[0014] The present invention aims to provide a 

method for analysing a backscatter signal for determining 
30 ±n an accurate way the spatial distribution of the beat 
length, the coupling length and the polarisation mode 
dispersion. 
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Summary of the invention 

[0015] The invention relates to a method for 

characterising an optical fibre link by its beat length, 
5 coupling length and polarisation mode dispersion 
distribution. It comprises the steps of 

- sending a pulsed signal along the optical fibre link and 
measuring the backscattered signal (preferably a POTDR 
signal) , after passing through a polariser, 

10 - deriving the length of the optical fibre, the average 
power difference between two successive minima of the 
backscattered signal and the number of maxima per unit 
length, 

- in an iterative way determining a beat length interval 
15 and an interval for the polarisation mode coupling 

parameter, until the length of the intervals is below a 
predetermined value, yielding a value for the beat 
length and the coupling length, 
calculating the polarisation mode dispersion. 
20 [001S] Preferably the POTDR signal is an ideal POTDR 

signa.1 . 

[0017] Advantageously the POTDR signal is the 

convolution of an ideal POTDR signal and a signal depending 
on th.e pulse shape. In an alternative embodiment the POTDR 
25 signal further is convoluted with a signal taking into 
account the time jitter introduced by the measurement set- 
up. 

[0018] In a specific embodiment a smoothing 

algorithm is applied to the POTDR signal. 
3 0 [0019] In another object the invention relates to a 

method for characterising an optical link consisting of a 
concatenation of several fibres, wherein the method as 
previously described is applied to each fibre. 
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[0020] The invention also relates to the use of a 

method as described above to locate the position of 
polarisation mode dispersion sources within an optical 
fibre link. 

5 [0021] In an advantageous embodiment the method as 

described is used in telecommunication networks . 
Alternatively, it can be applied in fibre sensing 
applications . 

10 Short description of the drawings 

[0022] Fig- 1 represents the set-up for measuring a 

POTDR trace. 

[0023] Fig. 2 represents three simulated POTDR 

traces for a given a and for different Lb values. 
15 [0024] Fig. 3 represents the number of maxima per 

unit length versus L B for different values of a. 
[0025] Fig. 4 represents simulated traces for a 

given beat length and different values of a. 

[0026] Fig. 5 represents the evolution of £ as a 

2 0 function of a for different values of L B . 

[0027] Fig. 6 represents the determination of the 

beat length interval . 

[0028] Fig. 7 represents the determination of the 

interval for a. 

2 5 [0029] Fig. 8 represents the distorsion of the POTDR 

signal due to the imperfections of the set-up (pulse shape 
+ time jitter) . 

[0030] Fig. 9 represents the evolution of n with the 

beat length when the imperfections of the set-up are taken 

3 0 into account . 

[0031] Fig. 10 represents the evolution of 4 with a 

when the imperfections of the set-up are taken into 
account . 
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[0032] Pig. 11 represents the simulated and measured 

probability density function of A. 

[0033] Pig- 12 represents the window algorithm that 

can be used for calculated the local PMD. 



Detailed description of the invention 

[0034] It is first shown how the polarisation mode 

dispersion from a POTDR signal can be determined in an 
accurate way (Lb is not approximately determined anymore) 

10 by analysing the statistical properties of its extrema . 
Suppose first that a falls within the range [0,15] 
degrees/m, which means that the coupling length is assumed 
to be larger than 12.8 m. This a interval has been chosen 
such that it corresponds to the values measured on real 

15 fibres . 

[0035] From a POTDR signal, the fibre length L, % 

and the number of maxima per unit of length, rim, are easily 
measurable. After measuring n m/ it is possible to determine 



a first interval \^B t > l b for the beat length of the 



20 fibre under test so that: 



L B e 



and are determined by the intersections of the 

curves n(L B/ o*) corresponding to a = 0 and a = 15 degrees/m 
with the straight line n = n^ as shown on Fig. 6. 
25 [0036] When the first beat length interval is 

determined, the curves ^(a f h B ) can be used to define an 

interval [^/ow>°up] for the a value of the fibre under test 
so that : 



ere 
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cr} ow and °)ip are determined by the intersections of the 



curves E(<7,Ii B ) for L B=^ B and L B"^ B with the straight 

*~ "low up 



line % = where £ ra is the £ value of the fibre under test 
measured from the POTDR signal. This principle is 
5 illustrated in Fig. 7. 

[0037] After the determination of o)ow and &up , the 

curves n(L B /Cr) for = &} 0 w and & = °\i P can be calculated. A 
new beat length interval is then deduced as previously 
explained by Fig. 6, but using the curves obtained for 

10 cr = cr} ow and cr = °\ip instead of a = 0 and a = 15 degrees/m. 

[0038] This new beat length interval denoted 



L 2 ,L\ 

B low B up 



\^ B low 9 ^ B up\ 



is included in -^r, and gives a smaller 



range for the beat length of the fibre under test. Fig. 7 
is then used again to calculate a new cr interval [ a Lj a «p] 
15 by using the curves 4(a,L B ) obtained for 

L B = l2 r and L B = L \ 

This new interval gives a smaller range for the a value of 
the fibre under test . This process is repeated until small 
enough intervals are obtained for the beat length L B and 

0 for a. 

[0039] When the beat length Lb and a are determined, 

the PMD is finally calculated by means of the relationship: 



PMD 



nr 1.545 



545 



{L B a 2 c 



<r*L 

<x 2 L 



0.439 



-1 + er 0439 



[0040] If the optical link to be characterised 

5 consists of a concatenation of several fibres, the method 



WO 2005/041449 



PCTYBE2004/000154 



11 

is applied to each fibre of the link and their beat and 
coupling lengths are therefore determined independently of 
each other- In this way f the PMD mapping is possible along 
the link. 

5 [0041] The method for analysing the backscatter 

signal as described offers several advantages : 

it requires only the access to one end of the optical 
fibre link to be characterised; 

it is based on a POTDR arrangement, but does not require 
10 the complete measurement of the complete state of 
polarisation of the backscattered signal . Only a simple 
polariser is required; 

it involves a POTDR measurement at only one wavelength; 
it gives accurate (not approximate) values for the beat 
15 length, the coupling length and the PMD. 

[0042] In practice, the final measured signal is 

also affected by noise. This noise is due to the detector 
noise and to the presence of a residual coherence noise. As 
the determination of the PMD is based on the detection of 
20 extrema, it is obvious that this noise will affect the 
results . 

[0043] As already mentioned before, the measured 

signal is also affected by the imperfections of the 
experimental set-up like the finite pulse shape and the 

25 time jitter of the electronic equipment. 

These phenomena (noise and imperfections of the measurement 
set-up) lead to a distortion of the POTDR signal. Therefore 
it cannot be used directly to deduce the PMD map. 
The width of the optical pulses launched into the fibre 

30 under test is finite and can be non-rectangular. 
Consequently, the backscattered power measured at the POTDR 
corresponds to an average of the ideal signal over a small 
spatial window depending on the pulse shape. The resulting 
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signal measured at the POTDR is in fact the convolution 
product between the ideal signal and the shape of the 
pulse : 

Ppulse(z) = P{z)®Pshape{z) 

5 where P pu ise(z) is the power measured by the POTDR according 
to a certain pulse shape. P(z) is the ideal POTDR signal 
and Pshape(z) is the power distribution along the pulse. 
Hence the Fourier transform of P pu ise(z) becomes: 

Ppuheik) = P{k)P shape (k) 

10 where P pu ise(k) , P(k) and P S ha Pe (k) are the Fourier transforms 
of P pu ise(z), P(z) and P S ha P e(z), respectively. 

[0044] Moreover the electronic equipment is not 

perfect: it is subject to a time jitter, i.e. the time 
delay between two successive pulses is not constant. This 
15 phenomenon strongly affects the averaging process of the 
POTDR. If we suppose that the number of pulses involved in 
the averaging process of the OTDR is large enough and that 
the time jitter is a random variable characterised by a 
Gaussian distribution of zero mean and standard deviation 
CFtj/ the signal detected at the POTDR can be written as: 



20 



z 2 



P m {z) = P pulse {z)®e 2a * 

where P ra (z) is the power measured by the POTDR taking into 
account the pulse shape and the time jitter. Hence the 
Fourier transform of P ra (z) becomes: 

jl 2 2 



25 Pmfr) = PpulseifyVtfe 2 

The pulse shape and the time jitter have for effect to 
distort the measured backscattered signal. An example is 
shown on Fig. 8, which has been obtained for parameters of 
the experimental set-up. The signals (a) and (b) correspond 
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to the ideal and affected signals, respectively. One can 
clearly observe that the signal is strongly distorted. 
[0045] Because of these imperfections, the method 

described above cannot be used directly: the curves n(L B ,a) 
5 and ^(a,L B ) are indeed based on ideal POTDR signals 
obtained by simulations. Therefore the effects of the pulse 
shape and the time jitter should be included in the 
calculation of n(L B ,a) and 2; (a, Lb) . In the case of the 
parameters of the experimental set-up and a pulse duration 

10 of 10 ns, these new curves are shown on Fig. 9 and 10. In 
order to deduce these curves, ideal POTDR traces have been 
simulated and the effects of the imperfections have been 
taken into account by using the above expressions for 
Ppuise(z) and P m (z) . Consequently, when the imperfections of 

15 the experimental set-up are taken into account, the PMD map 
of an optical link can still be determined. The same 
reasoning as the one explained previously has to be applied 
on the POTDR signal by using the new set of curves n(L B ,a) 
and £(a,Ls) described in fig. 9 and 10. As the pulse shape 

20 is now taken into account, it is possible to launch a 
larger pulse in the fibre (50 ns, for example), which will 
increase the system dynamics. The maximum measurable length 
will therefore be greater. 

[0046] As already indicated, the final measured 

25 signal is also affected by noise. This noise is due to the 
detector noise and the presence of a residual coherence 
noise. Its main effect is to add a series of minima and 
maxima to the POTDR trace. The values of n and £ cannot 
thus be correctly measured and the final calculated PMD 
30 will be erroneous. In order to solve this problem, a 
smoothing algorithm is applied to the backseat tered POTDR 
signal P m (z) . This algorithm consists of generating the 
signal P s (z) such that: 
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P s (z) is therefore the mean value of P m (z) on a distance 
interval L 9 around z. The number of maxima and the 
distances z k corresponding to the locations of the minima 
5 of the POTDR trace are then calculated from the smoothed 
signal in order to ignore the influence of the small peaks 
generated by the noise. Afterwards, when the positions z k 
are known, £ is calculated from P m (z) . When noise is 
present, the number of extrema of P s (z) depends of the 
10 smoothing parameter L s . So, the main problem is to choose 
the correct smoothing parameter to be applied to the POTDR 
signal . 

[0047] The location of the minima of the smoothed 

signal can be used to estimate the probability density 

15 function of A, denoting the difference in z between two 
successive minima. This estimation can then be compared to 
the density probability function of A of a POTDR trace 
obtained by means of simulations taking into account the 
parameters L B and a calculated from the smoothed signal. If 

20 the probability density functions of A obtained from P s (z) 
and from the simulation do not match, it means that the 
noise contribution is significant. The smoothing parameter 
L s will therefore be chosen such that the corresponding 
probability density function obtained from measurements can 

25 be correctly fitted by the probability density function 
derived from simulations: L s is therefore chosen such that 
the noise contribution is reduced. Fig. 11 illustrates this 
concept. Figures (a), (b) and (c) show the probability 
density functions obtained by measurement (circles) and by 

30 simulation (crosses) for a given optical fibre and for an 
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L s value of 3 , 4 and 5 m, respectively. One can clearly 
observe in this case L s = 3 m gives the best fitting. 
[0048] The invention also provides a technique which 

will allow to locate more accurately the position of the 
5 PMD source within the fibre. The basic idea consists in 
applying the described method on a distance window (1 km, 
for example) , which slides along the fibre length. By 
comparing the obtained PMD's on several successive windows 
as shown on Fig. 12, it is possible to locate more 

10 accurately high PMD sources within the optical fibre. 

[0049] POTDR can also be used in fibre sensing. The 

polarisation properties indeed depend on several external 
effects like strain, temperature, electric and magnetic 
fields. Strain, twist, temperature, electric and magnetic 

15 fields affect the birefringence within the fibre and 
consequently, the mean beat length. The method of the 
invention, which measures the spatial distribution of the 
beat length, therefore allows to measure the spatial 
distribution of these external effects after calibrating 

20 the sensing fibre. 



